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ABSTRACT 

We discuss the high enegry afterglow emission (including high energy photons, neutri- 
nos and cosmic rays) following the 2004 December 27 Giant Flare from SGR 1806-20. If 
the initial outflow is relativistic with a bulk Lorentz factor Tq ~ tens, the high-energy 
tail of the synchrotron emission from electrons in the forward shock region gives rise 
to a prominent sub-GeV emission, if the electron spectrum is hard enough and if the 
intial Lorentz factor is high enough. This signal could serve as a diagnosis of the initial 
Lorentz factor of the giant flare outflow. This component is potentially detectable by 
GLAST if a similar giant flare occurs in the GLAST era. With the available 10 MeV 
data, we constrain that Tq < 50 if the electron distribution is a single power law. For 
a broken power law distribution of electrons, a higher Tq is allowed. At energies higher 
than 1 GeV, the flux is lower because of a high energy cut off of the synchrotron emis- 
sion component. The synchrotron self-Compton emission component and the inverse 
Compton scattering component off the photons in the giant flare oscillation tail are 
also considered, but they are found not significant given a moderate Tq (e.g. < 10). 
The forward shock also accelerates cosmic rays to the maximum energy 10 eV, and 
generate neutrinos with a typical energy 10 14 eV through photomeson interaction with 
the X-ray tail photons. However, they are too weak to be detectable. 

Key words: satrs: neutron — stars: winds, outflows — hydrodynamics — Gamma 
Rays: bursts — accclcrtion of particles — elementary particles 



1 INTRODUCTION 

The Soft Gamma-ray Repeater (SGR) 1806-20 lies in the 
Galactic plane, at a distance of about Dl ~ 15.1kpc (Corbel 
& Eikenberry 2004; cf. Cameron et al. 2005). A giant flare 
originated from it on 2004 Dec. 27 is the brightest extra- 
solar transient event ever recorded (e.g., Hurley et al. 2005; 
Palmer et al. 2005). Radio follow-ups have resulted in de- 
tections of its afterglow (e.g., Cameron et al. 2005; Gaensler 
et al. 2005). Thanks to its brightness, an amazing variety 
of the data, including the source size, shape, polarization 
and flux at multi-frequencies as a function of time, have 
been collected (e.g., Gaensler et al. 2005; Cameron et al. 
2005; Gelfand et al. 2005). Even so, our understanding of 
the outflow is still in dispute. For example, the earliest af- 
terglow data obtained so far is about 7 days after the Giant 
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flare. At this epoch, even an initially relativistic outflow has 
been decelerated to the Newtonian phase by the interstellar 
medium (ISM). As a result, whether the outflow is relativis- 
tic initially (e.g. Wang et al. 2005; Dai et al. 2005) or not 
(e.g. Gelfand et al. 2005; Granot et al. 2005) is uncertain. In 
principle, similiar to the Gamma-ray Burst case (e.g., Krolik 
& Pier 1991), if the spectrum of the giant flare is nonther- 
mal, a lower bound of the initial Lorentz factor To ~ tens 
can be derived from the so-called "compactness argument" 
(e.g., Huang et al. 1998; Thompson & Duncan 2001; Nakar, 
Piran & Sari 2005; Ioka et al. 2005). Observationally, the 
giant flare spectrum may be thermal (Hurley et al. 2005) or 
nonthermal (Palmer et al. 2005), so that To could not be 
constrained well. 

In order to understand the dynamical evolution of the 
outflow better, early multi-wavelength (including optical 
and hard 7— ray band) observations are highly needed. The 
early optical emission has already been calculated in Cheng 
& Wang (2003) and Wang et al. (2005). In this work, we 
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focus on the high energy afterglow emission, including sub- 
GeV photons (see §2), high energy neutrinos and cosimic 
rays (see §3). High energy neutrinos from magnetars in the 
quiescent state have been discussed by Zhang et al. (2003) . 
Assuming the internal shock mechanism, the neutrino, cos- 
mic ray and TeV photon emission accompanying the prompt 
giant flare have been discussed recently (Ioka et al. 2005; 
Asano et al. 2005; Halzen et al. 2005). 



2 HIGH ENERGY PHOTON EMISSION 

We first take To = 10 as the typical Lorentz factor of the 
flow to do sample calculations. The effect of varying To will 
be discussed later. The isotropic energy of the outflow is 
taken as Ei BO ~ 10 46 ergs. In the following analytical discus- 
sion, we assume that the shocked electrons distributes as a 
single power-law dn/d^ c oc 7 t T p for 7 m < 7 C < 7m, where 
p ~ 2.5, 7m ~ 10 8 B' (B' is the shock generted mag- 
netic field strength, see equation (4)). Wang et al. (2005) 
find that a broken power-law distribution of electrons, i.e., 
dn/c?7e oc 7tT P1 for 7 m < j c < Jb and dn/d'fc oc 7,T P2 for 
7b < 7c < 7m, is required to interpret the chromatic radio 
afterglow lightcurve steepening around day 9. We therefore 
also include such a possibility in the numerical calculations 
(see §2.3). 

With the standard parameters, the relativistic outflow 
is decelerated by the ISM in a timescale 



t dcc « 300s Eh£ M ru 1/3 T 
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after which the ejecta moves with the Lorentz factor (for 

r>i/0j) 
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where n is the number density of the ISM, t obs is the observer 
time in unit of seconds. Throughout the work, we adopt the 
convention Q a = Q/10 a using cgs units. 

As usual, we assume e c and £b as the shock energy 
equipartition parameters for the shock accelerated electrons 
and the magnetic fields, repsectively. The minimum electron 
Lorentz factor reads 



7n 
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where C p = 3(p — 2) / (p— 1). The strength of shock generated 
magntic fields can be estimated as 



B' « 3.9 x 10~ 2 Gauss 4(- 2 ™o /2 [ r ( r - 1 )] 1/2 - 



(4) 



Throughout the work, the superscript ' represents the pa- 
rameter measured in the comoving frame of the ejecta. 

2.1 Inverse Compton Radiation 

A soft thermal X-ray tail emission modulated by the mag- 
netar period is typically detected after a giant flare hard 
spike. For the Dec. 27 event from SGR 1806-20, such a 
tail lasts for TUn ~ 300 with a typical photon energy 
ex ~ 30keV (e.g. Mazets et al. 2005) and a luminosity 
L x ~ 2 x 10 43 crgs s _1 (t obs /50s) _1 . For t ohs < T trai i, besides 
the synchrotron and synchrotron-self-Compton cooling pro- 
cesses (see 2.2 for detail), the electrons in the shocked region 
are also cooled by inverse Compton (IC) scattering off these 
X-ray tail photons. 



Since TUn is comparable to tdcc, the ejecta has not de- 
celerated significantly, i.e., T ~ To. In the comoving frame 
of the ejecta, the energy density of the X-ray tail reads 



0.27 ergs cm 3 Lx,^RisT 1 2 > 



(5) 
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where R is the radial distance of the forward-shock front 
from the central source. On the other hand, the magnetic 
energy density generated in the forward shock front reads 



Ub ~ 6 x 10 3 ergs cm 3 eB,-2r 2 no. 



(6) 



In the rest frame of the shocked electrons with a random 
Lorentz factor 7 e , the energy of the thermal tail 7 c ex/r is 
much larger than m e c 2 , so that the Klein-Nishina correction 
is important. For convience, we define x = 7 e ex/rm c c 2 ~ 
7e/17T. In the Klein-Nishina limit, aic = A(x)ot, where 

A( X ) EE f^f^-Kl + ^+iMl + to)-^], 

with the asymptotic limits A(x) ~ 1 — 2x + for 1 < 1, 
and A(x) w §a: _1 (ln2a; + \) for x » 1 (e.g. Rybicki & 
Lightman 1979). 

For illustration, we take t b s = Ttaii, at which R « 
2r 2 cT tai i « 1.8 x 10 15 cm r 2 T tai i, 2 . 5 . For 7e = 7m , we have 
x = 10.8, A(x = 10.8) « 0.1. The IC scattering is there- 
fore in the extreme Klein-Nishina limit, and the typical IC 
photon energy can be well approximated by 
1 ic 

hVrn 



hV^m e c « e e [(p- 2)/(p- l)](r - l)rm p c 2 
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The IC optical depth is 

r ~ A(W.8)a T nR/3 ~ 4 x 10" 



1 r> 

WO-H-15.26, 



(7) 
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so that the 10 GeV photon luminosity can be estimated by 

L10G0V ~ r(L x /e x )(/i^) 

~ 4.3 x 10 37 ergs s' 1 C p e c ,-o.5n -Ri5. 26^,2.5(9) 

where h is the Planck constant. For T ta n < tdcc, T ~ To, R oc 
tobs we have LiogcV oc Rt~ hs oc t® hB . We can then estimate 
the total number of the photons detectable by the Gamma- 
Ray Large Area Telescope (GLAST) 1 in construction 



iVtot(lOGeV) 



,A GL 



r 



£ioGcvdi bs]/10GeV 
0.03 T trail , 2 . 5 no(D L /15.1kpc)- 2 , (10) 



where Aglast ~ 8000cm 2 is the effective area of the 
GLAST. Since usually at least 5 photons are needed to 
claim a detection (e.g. Zhang & Meszaros 2001 and refer- 
ences therein), the above predicted JV to t is well below the 
threshold of GLAST. This component is undetectable for 
an energetic giant flare similar to the recent one even for a 
much closer SGR, for example, SGR 1900+14. 

The thermal tail photons would be also scattered by 
the electrons accelerated by the reverse shock. The reverse 
shock is expected sub-relativistic. At tdcc, 734 ~ 1.2, where 
734 is the Lorentz factor of shocked region relative to ini- 
tial unshocked outflow. Therefore, for the electrons accel- 
erated by the reverse shock, one has 7^ = e a [(p — 2)/(p — 
l)K m p/ m e)(734 — 1) ~ 37 by assuming the same parameters 
as in the forward shock region. Therefore e x will be scattered 
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to an energy ~ 7 m 2 e x ~ 30MeV. According to Eqs.(9) and 
(10), the detected number of photons essentially depends 
on the IC optical depth r and is independent on the typi- 
cal energy of the photons. We can then estimate the total 
number of the IC photons from the reverse shock region by 
comparing that in the forward shock region. First, the IC 
is now in the Thomson regime, i.e. aic — ot- Second, the 
total number of electrons contained in the reverse shock re- 
gion is about To times that in the forward shock region. The 
expected total number of the 30 MeV photons is therefore 

r„ 



jV to t(30MeV) 



-JV t ot(10GeV) ~ 3. 



(11) 



A(x = 10.8) 

The actual value should be smaller since the timescale of 
having a strong reverse shock could be shorter than Ttraii- 
Although this ~ 30MeV reverse shock component is more 
prominent than the ~ 10 GeV forward shock component, it 
is undetectable by GLAST, either. 

2.2 Synchrotron and Synchrotron-self-Compton 
Radiation 

For the forward shock emission, the cooling frequency i/ c , the 
typical synchrotron frequency v m and the maximum spectral 
flux F„ imax read (e.g. Cheng & Wang 2003; Wang et al. 2005) 

= 3.1 x 10 19 Hz Er^e-^noX^lil + Y)~ 2 , (12) 

(13) 



~< 2 P2 ,2 2 
y p- C/ iso,46 fc B,-2 fc c,-0.5'' o bs,3' 



F v , w = 474Jy S iso ,4 6 e^ 2 _ 2 nJ /2 ( T ^^)- 2 , (14) 

where Y is the inverse Compton parameter, which can be 
estimated by Y ~ [-1 + y/l + 4xe c /e B }/2 (e.g., Sari & Esin 
2001), where x = min{l, 2.67(7 m /7 c ) (p ~ 2) } is the radiation 
coefficient of the shocked electrons (see equation (A8) of 
Fan, Zhang & Wei (2005a)), and j c is the electron cooling 
Lorentz factor 

7.7 x 10 s 1 

7c ~ (1 + Y) TB'H ohs - (15) 

Notice that only the synchrotron self-Compton is consid- 
ered. The IC component discussed in §2.1 is in the extreme 
K-N regime at j c , giving a very small contribution to the Y 
parameter. So it is neglected. The resulting flux at a typical 
energy hu ohB — O.lGeV can be then estimated as 



P-2 P+2 
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k 0.1GeV 

For typical parameters, Y is in order of 1, so the flux is much 
higher than that of typical GRBs, so that this component 
can be detected by GLAST. Notice that there exists an upper 
limit on the synchrotron radiation energy 

hum ~ 2.8 x ufh^TB' ~ 1.2GeV V 1 , (17) 

above which a sharp cutoff is expected. 

Since the outflow is only mildly relativistic, the collima- 
tion effect is important in calculating the late lightcurves. 
Following Yamazaki et al. (2005), we adopt a half-openning 



0.3. This leads 
to a geometry corrected energy « _Ej so #?/4 w 2.2 x 
10 44 ergs E ,4e(^j/0.3) 2 , which matches that derived from 
the radio afterglow modeling (e.g., Wang et al. 2005). The 
above analytical calculations are only valid for T > 1/0 j. 
For T < 1/0 j , the jet sideways expansion effect is important. 
A rough estimate gives T(J B ) oc t~^ 2 , F„, max (J s ) oc t~^ s , 
f c (Js) oc const, and Vm(Js) oc f~ b 2 (e.g., Rhoads 1999), where 
J s represents the jet with important sideways expansion. 



/3 ,-,1/3 



-1/3 



Defining the "jet break" time t- s w 4500s r / E iso46 u 
for f bs > tj, the flux could be estimated as 
F Vobs ~ 1.8 x 10- 7 ergs cm" 2 GeV" 1 (t bsAj) _p (l + 
^"'(^obs/O.lGcV)- 15 / 2 . Since the typical Lorentz factor 
is small (~ 10), the above analytical treatment may not be 
a good approximation, and more detailed numerical calcula- 
tions are needed (see §2.3). The synchrotron-self-Component 
(SSC) luminosity (Lssc) could be estimated through the Y 
parameter, i.e. Y = Lssc/ L syn . This results in the maxi- 
mum SSC spectral flux (for idee < £obs < tj) 



F. 
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where is the typical SSC frequency 
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The resulting flux at fe^obs = O.lGeV reads 

^ob. = Kg™ (v obs /vl SC )- (p - 1)/2 

= 2.2 x 10 _8 ergs cm" 2 GeV -1 Y(l + y) p_3 Cp _1 
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For hv ohs < O.lGeV, this radiation component is much 
weaker than the synchrotron component. Beyond the syn- 
chrotron cutoff at hvohs ~ lGeV the SSC component domi- 
nates, but it is well below the GLAST threshold. 



2.3 Numerical results 

Similar to Huang et al. (2000) and Cheng & Wang (2003), 
we have calculated the dynamical evolution of the ejecta (see 
Fig.l) and the accompanying high energy photon emission 
(see Fig. 2) numerically. 

As shown in Fig. 1, the jet half-opening angle increases 
with time rapidly. With sideways expansion, the evolu- 
tion of the jet half openning angle could be written as 
(e.g. Huang et al. 2000) d8/dt obs = c s (F + W 2 - 1)/R, 
where c s w ^(4F + 3)(r 2 - l)/[3r(4r 2 - l)]c is the local 
sound speed. For T >• 1, this could be approximated as 
d6/dt h s ~ 1.2/(2rt bs). It is apparent that the sideways 
expansion of the jet is very important from the very begin- 
ning of the dynamical evolution if the initial Lorentz factor 
is as small as 10. As a result, there is no jet break in the 
(r— 1) lightcurve (Fig.l) the energy flux lightcurve (Fig. 2). 
This is different from the case of ultra-relativistic GRB out- 
flows, in which the sideways expansion is important only at 
later times. One conclusion drawn from Fig.l is that the 
ejecta accounting for the radio afterglow is nearly isotropic, 
which matches the observations well (e.g. Cameron et al. 
2005; Galensler et al. 2005). 
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According to Fig.2 where To = 10 is adopted, we can see 
that the predicted energy flux in the 0.05-0.15 GeV band is 
above the GLAST sensitivity (thick dashed line) , especially 
when a single power-law electron energy distribution (thin 
dashed line) is adopted. If the electron distribution is a bro- 
ken power law (solid line), the detectability by GLAST is 
only marginal. In the energy band above 1 GeV, the pre- 
dicted SSC energy flux (thin dash-dotted line) is always be- 
low the GLAST sensitivity (thick dash-dotted line), so that 
it is undetectable. 

In Fig. 3, we investigate the dependence of the predicted 
0.05-0.15 GeV energy flux (only the synchrotron radiation 
component is taken into account) on To. The general trend 
is that a higher To leads to a stronger sub-GeV emission. For 
To ~ tens, regardless of the distribution of the shocked elec- 
trons (single power law or broken power law), the predicted 
fluxes are all above the GLAST. For To ~ a few, only the 
single power law distribution model can yield to marginally 
observable 0.05-0.15 GeV photon emission. 

In principle, a measurement of the sub-GeV flux in 
the GLAST era could serve as a diagnosis of the intial 
Lorentz factor of the outflow. For SGR 1806-60, avail- 
able data already gives interesting constraints. Accord- 
ing to Matzets et al. (2005), the time averaged energy 
flux in the e 7 ~ 10 MeV band could be estimated as 
e 2 <£/V/de 7 ~ 1.6 x 10~ 6 ergs cm~ 2 s" 1 , where dN/de^ ~ 
10 _5 photons cm" 2 s" 1 keV" 1 is the photon number spec- 
trum of the tail emission at lOMeV. Comparing with our 
numerical results presented in Fig.3, the single power law 
distribution model with r = 50 is already above the ob- 
served level 2 . A Lorentz factor To > 50 is allowed only when 
a broken power law distribution of the electrons is assumed. 



3 COSMIC RAYS AND NEUTRINOS 

Below we estimate the maximum proton energy (ep 4 ) accel- 
erated by the forward shock. For simplicity, we only discuss 
W < tobs < tj. In general, e™ = min[e™(l), ef(2), e™(3)] 
satisfies three constraints (see also Fan, Zhang & Wei 
2005b). (1) The comoving shock acceleration time t a ~ 
tp/YeB' c should be smaller than the comoving wind ex- 
pansion time t' d ~ R/Tc, which yields e^(l) ~ eB'R. The 
numerical value reads 

M/-,\ r, , w inlf w 1/2 1/8^3/8 .-1/8 \ 

e p (1) ~ 2.4 x 10 eV e^_ 2 n ' Ej oA6 t ob l 3 . (21) 

(2) The comoving proton synchrotron cooling timescale 
t' cool — [67rmpC 3 /crTm 2 ]rep 1 B'~ 2 should be longer than the 
comoving acceleration timescale i a , which results in 

ef(2) * 2.6 x 10 21 eV ^E^^^ . (22) 

(3) The comoving proton cooling timescale due to photome- 
son interaction should also be longer than the comoving ac- 
celeration timescale t a . However, from equation (13), the 
typical frequency of the forward shock emission is too low 
to provide the target photons for photomeson interactions 

2 Our calculated energy flux is in the e 7 = 50 — 150 MeV band. 

However, since ^obs^obs ^ u ^ob s P '^ 2 verv weakly depends on 
v obs> the results could approximately apply to the 10 MeV band 
as well. 
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Figure 1. Dynamical evolution of the ejecta. The dotted line 
represents T — 1; the solid line represents (multiplied by 10); 
the dash-dotted line represents the SSC parameter Y (multi- 
plied by 10); and the dashed line represents the cut off en- 
ergy of synchrotron radiation hvyy (normalized to 0.12GeV), all 
as functions of time. Following initial parameters are adopted. 
E iso = 10 46 ergs, 6, = 0.3, T = 10, n = lem- 3 , e c = 0.3, 
ee = 0.01 and p = 2.5. The starting point of the calculation 
is taken as Rq = 10 13 cm. 



at the A resonance, so the effect of photomeson interaction 
is too small to change the proton cooling process. 

Therefore, one has ef = ef(l) ~ 2.4 x 10 17 eV. The 
source location of SGR 1806-20 is about 10° from the Galac- 
tic center. In the region near the Galactic center the mag- 
netic field structure is poorly constrained. The time delay 
due to the interstellar random magnetic field can be approx- 
imated as (e.g. Asano et al. 2005) 

T dclay ~ (eB G D L / eC R) 2 (//c), (23) 

where Bg ~ 10~ 6 G is the average magnetic field strength 
in the Galaxy, ecR is the typical cosmic ray energy, and 
Z ~ 10 — lOOpc is the correlation length of the magnetic 
field (e.g., Asano et al. 2005). One then gets T do i ay r; 6 x 
10 5 yr B 2 }j _ 6 (D L /15.1kpc) 2 (//10pc)e c : R , il 7- As a result, these 
cosmic rays become a part of the cosmic ray background. 

As shown in equation (13), the typical frequency of the 
forwardshock emission is too low to provide the target pho- 
tons for the photomeson interaction at the A-resonance. The 
only interesting source of the neutrino emission is then the 
photomeson interaction during the early epoch when the X- 
ray tail overlaps with the shocked region. In the comoving 
frame of the ejecta, the thermal tail photons with energy 
w ex/r interact with the protons with energy 

e p ~ 0.3r 2 GeV 2 /e x ~ 10 16 eV rf (ex/30keV) _1 . (24) 
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Figure 2. Predicted photon energy fluxes as compared with the 
GLAST flux sensitivity. The thin dashed line and the solid line 
represent the flux in the energy range 0.05-0.15 GeV for differ- 
ent electron energy distributions, where "S" denotes single power 
law and "B" denotes broken power law. Only synchrotron com- 
ponent is calculated since the SSC component is much dimmer. 
The thin dash-dotted line represents the flux in the 1-300 GeV 
energy band. Only the SSC component is calculated since this is 
above the synchrotron cutoff energy. The thick dashed line and 
the thick dash-dotted line represent the GLAST threshold in the 
energy range 0.05-0.15 GeV and 1-300 GeV, respectively. The 
GLAST threshold is defined by requiring that during the inte- 
gration timescale ~ t b s at least 5 photons arc collected. For the 
thin dashed and dash-dotted lines, same parameters as those in 
Fig.l are taken. For the solid line, the parameters are the same as 
those taken in Fig.l except that pi = 2.2, p2 = 3.2 (rather than 
p), and 7), /7m = 120 are adopted. 

These protons lose ~ 20% of their energy at each jry inter- 
action, dominated by the A-resonance. Approximately, half 
of pions are charged and decay into high energy neutrinos 
7t + — > fi + + — > e + + v c + + Vp, with the energy dis- 
tributed roughly equally among the decay products (e.g., 
Ioka et al. 2005). Therefore the neutrino energy is ~ 5% of 
the proton energy, i.e., 

El ,~5x 10 14 eV r 2 (e x /30keV) _1 . (25) 

The comoving number density of the thermal photons at the 
radius R ~ 10 15 cm is 

nx w E&/(ex/r) m 5.5xl0 7 ix,4 3 i?i 5 2 rr 1 (ex/30kcV)- 1 .(26) 

The fraction of the energy converted to pions can be es- 
timated by the number of the p — 7 interactions occuring 
within the shock with the characteristic width AT? ~ R/T, 
i.e. 

f n ~ 0.2n x a A R/T 

~ 5.5 x lO^Lx^r/rr^ex/SOkeVr 1 . (27) 



Figure 3. The dependence of predicted 0.05-0.15 GeV photon 
energy flux on the initial Lorcntz factor To, whose actual values 
have been marked in the figure. Both single power law ("S") and 
broken power law ("B") distributions of the electrons are calcu- 
lated. The thick dashed line represents the GLAST sensitivity in 
the energy range 0.05-0.15 GeV. Except To, other parameters are 
the same as those taken in Fig. 2. 

where <ta ~ 5 x 10 _28 cm 2 is the cross section of the 
A— resonance. For a neutrino detector with an area A dct ~ 
10 10 cm 2 , the expected event number is 

N v ~ P„^UA det E iso /(32nDle„) ~7x 10" 5 , (28) 

where P„^ M ~ 3.5 x 10~ 4 (e„/10 15 cV) ' 5 is the probability 
that a neutrino produces a detectable high energy muon 
for e„ > 10 3 TeV. We can see that the predicted neutrino 
number is well below the detection threshold of the most 
powerful neutrino detectors under construction. The main 
reason is that compared with GRBs, f w (Eq.[27]) is much 
smaller. 



4 SUMMARY 

We show that if a giant flare similar to the 2004 Dec. 27 event 
happens in the GLAST era, a strong sub-GeV flare shortly 
after the flare (originated from the hard tail of the syn- 
chrotron emission from the forward shock region) should be 
detectable if the outflow is relativistic. A positive/negative 
detection in the sub-GeV band would then give a diagnosis 
of the initial Lorentz factor of the outflow. With the avail- 
able 10 MeV data (Mazets et al. 2005), we constrain T < 50 
for the Dec. 27 event if the electron distribution is a single 
power law, although a higher To is allowed if the electron 
distribution is a broken power law. At higher energies (e.g. 
above 1 GeV), a cutoff of the synchrotron emission is ex- 
pected. Neither the synchrotron self-Compton emission in 
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the forward shock region, nor the inverse Compton off the 
X-ray tail emission, could give a detectable flux for GLAST. 

The forward shock is able to accelerate protons to an 
energy ~ 10 17 eV. But the time delay for these cosmic rays 
to reach us is very long, i.e. ~ 10 6 years. Neutrinos with 
an energy 10 14 eV are also predicted, but the flux is too low 
to be detected. Therefore, for a giant flare similar to the 
Dec. 27 event taking place in the GLAST era, the most, and 
perhaps the only, interesting high energy afterglow emission 
is the bright sub-GeV photon emission lasting for thousands 
of seconds. 
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